Nonmonotonical crossover of the effective susceptibility exponent 
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We have numerically determined the behavior of the magnetic susceptibility upon approach of the 
critical point in two-dimensional spin systems with an interaction range that was varied over nearly 
two orders of magnitude. The full crossover from classical to Ising-like critical behavior, spanning 
several decades in the reduced temperature, could be observed. Our results convincingly show that 
the effective susceptibility exponent 7cff changes nonmonotonically from its classical to its Ising value 
when approaching the critical point in the ordered phase. In the disordered phase the behavior is 
monotonic. Furthermore the hypothesis that the crossover function is universal is supported. 
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At a continuous phase transition several thermody- 
namic observables diverge as a power of the temperature 
distance to the critical point. These powers, or critical 
exponents, have universal values which are identical for 
large classes of systems. For example, uniaxial ferromag- 
nets, binary alloys, simple fluids, binary mixtures, ionic 
solutions, and polymer mixtures all belong to the three- 
dimensional Ising universality class. However, the cor- 
responding power-law behavior is only observed asymp- 
totically close to the critical point. As stated by the 
Ginzburg criterion [Q, classical or mean-field-like criti- 
cal behavior may be observed at temperatures further 
away from the critical temperature T^. The explanation 
of this crossover in terms of competing flxed points of a 
renormalization-group transformation is one of the great 
achievements of Wilson's renormalization theory. Never- 
theless, the precise nature of the crossover between these 
two universality classes is still subject to debate. The- 
oretically several attempts have been made to approxi- 
mately calculate crossover functions. For instance, NicoU 
and Bhattacharjee |Q solved the renormalization equa- 
tions in d dimensions to second order in e = 4 — by 
applying a specific matching condition, whereas Bagnuls 
and Bervillier used massive field theory in d = 3. The 
results of Belyakov and Kiselev Q are phenomenological 
generalizations of first-order e-expansions. All these re- 
sults are only valid in the symmetric phase (T > Tc) and 
suggest that the crossover behavior is universal. How- 
ever, Anisimov et al. Q claimed that, while at criticality 
microscopic cutoff effects may be neglected compared to 
the infinite correlation length, this is no longer the case 
in the crossover region. This implies that the crossover 
functions cannot be represented as universal functions of 
one variable. A particular question concerns the vari- 
ation of the so-called effective exponents describing the 
continuous change from one type of power-law behavior 
to another in the crossover region. Whereas some cal- 
culations predict a strictly monotonical variation, others 
indicate that a nonmonotonical variation might be possi- 



ble. While on the theoretical side several important open 
questions remain, the experimental situation is hardly 
better. Measurements in the critical region are difficult 
and accurate results are scarce. Fisher Q has discussed 
experiments on micellar solutions (expected to belong to 
the Ising universality class) [Q that yielded values for the 
susceptibility exponent 7 that lie below the classical value 
7mf = 1, while the Ising value is given by 71 = 1.237 [|. 
He argues that these results can be incorporated in a 
standard scaling description of crossover behavior if one 
allows for an effective susceptibility exponent that varies 
nonmonotonically as a function of the reduced tempera- 
ture t = [T — Tc)/Tc. More recently, Anisimov et al. |^ 
measured a susceptibility x for which they found that the 
logarithmic derivative 7off = — dlnx/rfln|t| approached 
the Ising value from above upon approach of the critical 
point. As this implies a nonmonotonical variation of joS, 
Bagnuls and Bervillier subsequently suggested that the 
measurements might have been taken outside the critical 
region, see Refs. pO| , pT| . Indeed, since the crossover re- 
gion is expected to span several decades in the reduced 
temperature |^,^ , in many experiments the full crossover 
behavior cannot be observed. At the same time, the large 
extent of the crossover region reinforces its experimental 
relevance: many measurements of critical exponents are 
actually made within the crossover region and thus only a 
detailed knowledge of the crossover behavior guarantees 
a correct interpretation of the data. 

Although it is tempting to apply numerical simulations 
to shed some light on these issues, in practice one en- 
counters difficulties comparable to those experienced by 
experimentalists. In particular the size of the crossover 
region constitutes a towering hurdle. A major effort has 
been undertaken in Ref. jl^ for three-dimensional poly- 
mer mixtures, in which crossover occurs as a function of 
the polymer chain length. These systems offer the advan- 
tage that the crossover can be influenced both by vary- 
ing the temperature and by changing the chain length. 
Despite chain lengths of up to 512 monomers, the re- 
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suits did not span the full crossover region. Mon and 
Binder [ p^ examined the two-dimensional Ising model 
with an extended range of interaction R, where crossover 
from Ising to classical behavior occurs when R is in- 
creased. They studied crossover in finite systems at 
T — Tc- Even in these systems the mean- field regime 
could hardly be reached. In Ref. we showed that a 
new Monte Carlo (MC) cluster algorithm for long-range 
interactions could be applied to this model, leading 
to a speed increase of many orders of magnitude com- 
pared to conventional algorithms. 

In this Letter we use this algorithm to study two- 
dimensional Ising systems with a variable interaction 
range and present results for the crossover behavior of 
the magnetic susceptibility at temperatures below and 
above Tc. Although two-dimensional systems are simpler 
than their three-dimensional counterparts, this model ex- 
hibits a surprising behavior. In particular, a qualitative 
difference between T < and T > Tc is found. The ad- 
vantage of examining two-dimensional instead of three- 
dimensional systems is the much larger variation of the 
critical exponents in the crossover region and the accessi- 
bility of larger interaction ranges, which makes it feasible 
to cover the full crossover region. 

The model under investigation was introduced in 
Ref. and is defined by the following Hamiltonian, 

n/kBT=-J2Kd{rr-rj)s,.3, , (1) 

where the spins s take the values ±1, the sum runs over 
all spin pairs, and the spin-spin coupling depends on 
the distance |r| between the spins as Kd{r) — cR:^'^ for 
|r| < Rm and Kdir) = for |r| > R^- For finite Rm 
the critical behavior of this model will be Ising-like, but 
for Rm — > oo it will be classical. This implies a singular 
dependence of the critical amplitudes on i?^, which was 
first derived on phenomenological grounds in Ref. [ p^ . 
In Ref. ijlj] a renormalization derivation of these singu- 
lar dependences was given, which in addition revealed 
logarithmic corrections for rf = 2. To avoid lattice ef- 
fects we formulate range dependences in terms of an 
effective interaction range R, which is directly related 
to Rm The Ginzburg criterion introduces a param- 
eter G oc i?~2d/(4-d) j^^j^g Ginzburg number) which de- 
termines whether the critical behavior will be Ising-like 
(t ^ G) or classical {t » G). In the latter case, care 
must be taken that t is still within the critical region. For 
many experimental systems the Ginzburg number is not 
small and one has left the critical region before observ- 
ing the full crossover to classical critical behavior. In our 
model system, G is adjustable so that we can vary t/G 
over the full crossover region while keeping t sufficiently 
small. On the other hand, for a too small Ginzburg num- 
ber the critical point must be approached very closely to 
access the Ising regime. The diverging correlation length 



is then in our simulations truncated by the finite system 
size L. Therefore we construct the crossover function by 
studying systems with various values of G (interaction 
ranges) such that t has to be varied only within a lim- 
ited range (but in such a way that the results for several 
different G overlap at fixed t/G). 

We have carried out MC simulations of square sys- 
tems with periodic boundary conditions containing up to 
1000 X 1000 spins in which each spin interacts with up 
to 31416 neighbors. This corresponds to an effective in- 
teraction range R of 71 lattice spacings or intermolecular 
distances. To avoid systematic errors in the determi- 
nation of the crossover behavior, an accurate estimate 
of Tc as a function of R is required. For systems with 
interaction ranges up to i? « 8.3 results for Tc(i?) have 
been obtained in Ref. ||l^. For larger ranges the critical 
temperature can be calculated to a comparable accuracy 
from a renormalization expression for Tc{R) . Further 
simulational details will be presented elsewhere |p^ . 

In the two-dimensional Ising model 71 — 7/4 and 
the susceptibility x diverges for t f as (— t)^*"/* 
and for t | as A^t~''/^. The critical amplitudes 
for the nearest-neighbor model are known exactly p7| , 
= 0.025537 ... and A+ ^ 0.96258 .... Note the very 
large asymmetry, Af /A^ ^ 38. Mean-field theory pre- 
dicts a susceptibility that for f t diverges as l/{—2t) 
and for i J, as 1/t, i.e. a susceptibility exponent 7mf = 1 
and a much smaller ratio A^-p/A^-p — 2. As derived 
in Refs. |]l^,Q, the Ising critical amplitude of the sus- 
ceptibility is proportional to R^^/^. Thus, in a graph 
displaying the results for various ranges as a function 
of the crossover variable t/G cx tR^ a data collapse is 
obtained for x/R?. The susceptibility is related to the 
average magnetization per spin m. In our simulations we 
have for i < sampled the connected susceptibility given 
by the fluctuation relation x — L'^({m^) — {\m\)'^) / k-gT , 
whereas for i > we have used x = L'^ {rri^) / k-gT . 

In Fig. |l| we show the magnetic susceptibility below Tc 
for various system sizes and interaction ranges. This 
graph exhibits several interesting features. For very small 
values of |t| the curves lie almost horizontal; this is the 
finite-size regime where the correlation length is trun- 
cated by the system size. For somewhat larger values 
of |t| the curves start following the Ising asymptote with 
slope (i.e., the logarithmic derivative dlnx/dln \t\) —7/4. 
This is the critical behavior as it is experimentally mea- 
sured close to Tc. Also the critical amplitude A^ is accu- 
rately reproduced by the simulations. At even lower T, 
we see that the curves for systems with small interaction 
ranges start to deviate from the Ising asymptote toward 
the mean- field asymptote (slope —1) without actually 
reaching it. These systems have left the critical region 
and the order parameter shows strong saturation effects 
(which decreases the susceptibility). However, systems 
with a larger interaction range clearly cross over to the 
mean-field asymptote, reproducing the mean-field criti- 
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FIG. 1. Crossover behavior of the connected susceptibil- 
ity X for various ranges and system sizes. In this and all 
following figures the numbers in the key refer to values for 
the interaction range R. 



cal amplitude ^^f- ^'^^ even lower temperatures these 
systems also exhibit saturation effects, which for R > 8.3 
are accurately described by mean-field theory (dashed 
curves). However, the outstanding feature of this graph 
is the region between the Ising and the mean-field asymp- 
tote. Namely, before settling at the latter asymptote, 
the curve describing the susceptibility first has (in this 
double-logarithmic plot) a slope that is less steep than 
in the mean-field regime. That is, 7^ < 1 (the super- 
script minus sign indicates that we are considering the 
case t < 0). To illustrate this effect more clearly we have 
reproduced Fig. ^ without the data that are plagued by 
finite-size effects or lie outside the critical region. Fur- 
thermore, we have corrected for the saturation effects for 
systems with large ranges. While this is a real physical 
effect, it can be removed by applying a correction fac- 
tor accounting for the difference between the asymptotic 
mean-field susceptibility and the mean-field susceptibility 
affected by saturation. The resulting graph is shown in 
Fig. 0. The nonmonotonical variation of the slope is now 



clearly visible. The data for different interaction ranges 
1 ^ ^ ^ 70 overlap for considerable intervals of tR^ . The 
perfect collapse of these data lends strong support to the 
hypothesis that the crossover curve is universal and spans 
several decades in the reduced temperature. In addition, 
it follows from the renormalization treatment in Ref. |l^] 
that the correlation length ^ decreases as t~'^ with an 
amplitude which is for d = 2 to leading order indepen- 
dent of R. Thus, at a fixed value of tR^ the curves for 
different ranges have different values for ^ and the fact 
that they collapse implies that the ratio between ^ and 
the lattice spacing a does not affect the crossover curve. 
This is markedly different from the results of Ref. § for 
d = 3. Also the influence of irrelevant fields is not visible 
in the data collapse. To connect to experimental results, 
we have plotted the effective exponent 7^ (obtained by 
numerical differentiation) in Fig. ^ Starting from Tc, 7^ 
first steeply decreases to a minimum below 7mf and then 
gradually rises to the asymptotic mean-field value. 

Now we turn to the symmetric phase, for which a data 
collapse of the susceptibility x is shown in Fig. ||. Just as 
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FIG. 3. The effective susceptibility exponent 7^,^ below Tc, 
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FIG. 2. Crossover curve for the connected susceptibility x- 
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FIG. 4. Crossover behavior of the susceptibility x for vari- 
ous ranges and system sizes. 
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below Tc, finite-size effects occur for very small t. Outside 
the finite-size regime the data for various R nicely col- 
lapse on the Ising asymptote, again with slope —7/4 but 
with a much larger amplitude. For higher temperatures 
the curves appear to gradually approach the mean-field 
asymptote. However, only for larger interaction ranges 
the critical amplitude is reproduced. This stresses 
an important point: above no saturation of the or- 
der parameter occurs, marking the end of the critical 
region, but the system smoothly passes over to regular 
(noncritical) behavior. In this high-temperature region 
the susceptibility decreases proportional to l/T. It is 
this behavior that is seen in the graph at high tempera- 
tures for systems with small interaction ranges. As right- 
fully stressed by Bagnuls and Bervillier [|o|, this behav- 
ior should be clearly distinguished from classical critical 
behavior. Thus, it is by no means disturbing that the 
curves for small R deviate from the mean-field asymp- 
tote and this does not imply a nonuniversal character of 
the crossover curve in the critical region. Disregarding 
these systems that have left the critical region, we note 
that above Tc the susceptibility smoothly crosses over 
from Ising-like to classical critical behavior and that the 
effective exponent 7^ decreases monotonically from 7/4 
toward 1, as visualized in Fig. 0. 
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FIG. 5. The effective susceptibility exponent 7^ above Tc. 

In conclusion, we have presented crossover curves 
for the magnetic susceptibility in two-dimensional Ising 
models with medium-range interactions, both for i < 
and i > 0. Unlike previous treatments, which all suffered 
from some systematic limitations (like extrapolating low- 
order £-expansions to physical dimensions), the present 
approach for the first time gives an explicit description 
of crossover scaling functions for critical phenomena. At 
least in principle, these curves could be directly compared 
to experimental results. The large interaction ranges that 
could be accessed allowed us to observe the crossover be- 
tween Ising-like and classical critical behavior for nearly 
six decades in the crossover variable. This has yielded 
for the first time strong numerical evidence that below Tc 



the effective susceptibility exponent 7efi varies nonmono- 
tonically between its Ising value and its classical value. 
Above Tc, on the other hand, the exponent shows a mono- 
tonical variation between the two limiting values. Thus, 
the plausibility of the occurrence of a minimum in 7off in 
three-dimensional systems, at least in the phase of broken 
symmetry, has been greatly increased. Furthermore, the 
fact that the crossover curves for many different inter- 
action ranges collapse supports the hypothesis that this 
curve is universal. 
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